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Abstract 

The magnetic field dependent critical current density j c (B) of a MgB2 bulk 
sample has been obtained by means of magnetization hysteresis measure- 
ments. The jc(B) curves at different temperatures demonstrate a crossover 
from single vortex pinning to small-bundle vortex pinning, when the field 
is larger than the crossover field B s ^. The temperature dependence of the 
crossover field B s ]j(T) is in agreement with a model of randomly distributed 
weak pinning centers via the spatial fluctuations of the transition temper- 
ature (<5T c -pinning) , while pinning due to the mean free path fluctuations 
((^-pinning) is not observed. 
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The recent discovery of superconductivity in the intermetallic compound MgB2 [[jj with 
transition temperature at 39K has led to intensive experimental and theoretical activities 
P^T5|, with the purpose of understanding the basic mechanism of superconductivity and 



the vortex pinning mechanism governing the critical current density j c in this new supercon- 
ductor. Although the critical current density has been improved greatly since its discovery, 
the underlying pinning mechanism is still under investigation. 

In type-II superconductors, the most important elementary interactions between vor- 
tices and pinning centers are the magnetic interaction and the core interaction |16|^2|. The 
magnetic interaction arises from the interaction of surfaces between superconducting and 
non-superconducting material parallel to the applied magnetic field. In technical type-II 
superconductors with a high Ginzburg-Landau (GL) parameter k, the magnetic interaction 
is usually very small and disappears with increasing magnetic field. The core interaction 
is usually more effective in technical type-II superconductors due to the short coherence 
length and the larger penetration depth (high k). This interaction arises from the coupling 
of the locally distorted superconducting properties with the periodic variation of the super- 
conducting order parameter. Two mechanisms of core pinning are predominant in type-II 
superconductors, i.e. 5T C - and ^/-pinning. Whereas <5T c -pinning is caused by the spatial 
variation of the GL coefficient a associated with disorder in the transition temperature T c , 
variations in the charge carrier mean free path I near lattice defects are the main cause of 
^/-pinning. 



It has been reported by Griessen et al. (T^j that the ^/-pinning mechanism is dominant 



in both YBa2Cu307 and YE^C^Og thin films. For the new superconductor MgB 2 , a high 
k value of 26 has been reported flO |, it is therefore expected that the magnetic interaction 



is negligible, while the core interaction is more important. However, it has not been experi- 
mentally determined whether the ^/-pinning or the 5T c -pinning is the dominant mechanism 
in MgB2- The purpose of this Letter is to report measurements of the critical current den- 
sity of this new material to achieve an understanding of the vortex pinning mechanism and 
to demonstrate that in MgB2 governed by bulk pinning, 6T C is the only important pinning 



mechanism. 

All measurements have been performed on a MgB 2 bulk sample, which was prepared by 
conventional solid state reaction |2"3]| . High purity Mg and B (amorphous) with a nominal 
composition ratio of Mg:B=1.2:2 were mixed and finely ground, then pressed into pellets 10 
mm in diameter with 1-2 mm thickness. Extra Mg was added in order to make up for loss of 
Mg at high temperatures. These pellets were placed on an iron plate and covered with iron 
foil, then put into a tube furnace. The samples were sintered at temperatures between 700 
and 1000°C for l-14h. A high purity Ar gas flow was maintained throughout the sintering 
process. A sample with T c = 38.6 K and dimensions of 2.18 x 2.76 x 1.88 mm 3 was cut from 



the pellet. Phase purity was determined by XRD ]24| and grain size by SEM. Only a small 
level of MgO (less than 10%) was found and the grain size was determined to be about 200 
/iin. 

Fig. |l| shows the magnetization hysteresis loops of the MgB 2 sample every 2 K in the 

are 



14-36 K range. The results at lower temperatures, which have large flux jumping |23 
not shown here. The symmetric magnetization hysteresis loops with respect to the magnetic 
field indicate the dominance of bulk current up to temperatures near T c , rather than surface 
shielding current. Therefore, the bulk pinning is dominant in this sample, while the surface 
pinning is negligible. As a comparison, we show in the inset of Fig. [l] the magnetization 
hysteresis loop of a pressed MgB 2 sample at 5 K. This sample is fabricated by pressing the 
MgB 2 powder into a pellet without sintering. The loop is highly asymmetric, showing a 
large reversible magnetization resulting from the surface current, and surface pinning plays 
an important role in this sample. No flux jumping is observed down to T=5 K, indicating 
that the grains in the sample are decoupled. The surface pinning effect has also been observed 
by Takano et al. [13| in their powder sample and bulk sample sintered at low temperature. 

From these M(H) loops, we can calculate the critical current density using j c = AM/a(l — 
a/36), with a, b the width and length of the sample perpendicular to the applied magnetic 
field, respectively. The resultant j c (B) curves at various temperatures are shown in a double 
logarithmic plot in Fig. as different symbols. As can be seen from the plateau at low 



magnetic field, j c initially has a weak dependence on the field. When the magnetic field is 
increased beyond a crossover field, it then begins to decrease quickly. The crossover field 
decreases with increasing temperature. Further increasing the magnetic field results in a 
faster drop in j c near the irreversibility line, which is obtained by using a criterion for the 
critical current density j c = 100 A/cm 2 . The results are shown as open circles in Fig. |3|. 
The best fitting of the data yields the result, 



B irr (t) = B iTT (0) (1 - t 2 f , (1) 

shown as solid line in Fig. |3|, with t = T/T c . For high temperature superconductors, a 
(1 — t) 3 / 2 behavior is usually observed |25| and has been explained by means of giant flux 
creep p6|j27[| . The (1 — t) 3 ^ 2 law is actually an approximation of the (1 — 1 2 ) 3//2 law as t — > 1, 
because (1-t 2 ) « (l-t)(l + t) « 2(1 — t). Also plotted in the fi gure is the B c2 (T) data taken 
from Takano's work (the dashed line is just guide to the eye), showing that B irv is well below 
B c2 . Therefore, giant flux creep also plays an important role in the new superconductor 
MgB 2 . 

j c {B) characteristics very similar to those shown in Fig. ^] have been observed by other 
groups P JT^JT3[ ]. Based on different physical assumptions on summation of the elementary 
pinning force f p to obtain the macroscopic pinning force F p , different pinning models yield 
different j c {B) characteristics. The simplest model is the direct summation of f p to have 
F p = j c B = n p f p , where n p is the density of pinning centers in the sample. This strong 
pinning model neglects the influence of the flux line lattice. When the influence is taken into 
account, one have F p = j c B = n p f 2 (u / f p )d/al, where u is the maximum distortion of the 
flux line lattice, d is the range of the pinning force, typically of the order of xi, and ao the flux 
line lattice constant. These two strong models yield j c {B) characteristics of j c oc B^ 1 and 
j c oc 5~ a5 , respectively. Due to the large densities of the pins n p and the small elementary 
interaction forces f p , these two models are not representative for most real pinning systems 
|28fl. For randomly distributed weak pinning centers, the macroscopic pinning force F p can 



be estimated using the basic concept of collective pinning |29| , which has been proved to be 



very successful in most real pinning systems, 



with the correlation volume V c = R 2 C L C) with the correlation lengths R c perpendicular to 
the field direction and L c along the vortex line, and the pinning parameter W = n p < f£ >. 
R c and L c depend on the applied magnetic field, the dimension of the flux line lattice (3D 
or 2D), and the elasticity or plasticity of the flux line lattice. For the 3D elastic flux line 



lattice, it has been derived by Blatter et al. |17| that j c is field-independent when the applied 



magnetic field is lower than the crossover field B s ^ (single vortex pinning) 

B s b = Psb—B c2 (3) 
Jo 

where /3 s b ~ 5 is a constant, j = 4B c /3^/6fj>o\ the depairing current, B c = $ /2v^2VA£ the 
thermodynamic critical field, B c2 = /io$o/2vr^ 2 the upper critical field and j sv the critical 
current density in the single vortex pinning regime. When the applied field is larger than 
-B s b (small bundle pinning), j c (B) follows a exponential law 



j c (B) « j c (0) exp 



U 



(4) 



When the applied magnetic field is larger than B\ h = Ab-B C 2(jsv/jo)[hi(« 2 j S v/jo)] 2//3 (where 
/?ib is a constant w 2), this large bundle pinning regime is governed by a power law j c {B) oc 
B~ 3 . The 2D elastic flux line lattice shows single pancake pinning at low magnetic fields 
with field-independent j c , then a 2D collective-pinning region at higher fields with j c oc 1/B. 
For fields higher than a crossover field B 3D , three-dimensional pinning is predicted. 

From the characteristics of the j c {B) curves shown in Fig. |2|, it is expected that the 
3D elastic pinning model (single vortex pinning followed by small bundle pinning, then 
large bundle pinning) may the dominant pinning mechanism in MgB2. We therefore use 
Eq.(f|) to fit the j c {B) curves, with fitting parameters j c (0) and B . The fitting results for 
different temperatures are shown as solid lines in Fig. |2|. At intermediate fields, Eq.([|) fits 
the experimental data very well, while deviations from the fitting curves can be observed at 



both low and high fields. A clearer plot is shown in Fig. |4], where the j c {B) curve at 24 K 



is shown in a double-logarithmic plot of —log\j/j(B = 0)] versus the applied field, which 
clearly shows a straight line at intermediate magnetic fields. The deviation at low fields is 
denoted as _B s b, indicating the crossover from the single vortex pinning regime to the small 
bundle pinning regime. The point of deviation at high fields was first considered as the 
crossover field from small bundle pinning to large bundle pinning. However, when we fit the 
j c (B) data at high fields to the power law j c {B) oc B~ n , n is found to be as large as 20 rather 
than the theoretically predicted value of 3, making it unlikely that the system changes to the 
large bundle pinning regime. As the high field deviation is very close to the irreversibility 
line, which results from giant flux creep, it is likely that the deviation at high field may 
result from large thermal fluctuations, which lead to the rapid decrease in j c , and therefore 
is denoted as _B t h> indicating thermal fluctuations. Other j c {B) expressions have also been 
tested for the fitting, such as j c (B) oc 1/[1 + (B/B ) n } and j c (B) oc exp[—(B/B )], but both 
yield poor fitting results. The inset of Fig. ^ shows the j c {B) of the pressed sample at 5 K 
in a double- logarithmic plot, which indicates that when the surface pinning is important, 
the exponential drop in j c (B) [see Eq.(|])] no longer applies, but a power law j c {B) oc B~ 1,2 



The crossover field B& we have obtained between single vortex pinning and small bundle 
pinning as a function of temperature is shown in Fig. |5| as open circles. We now compare 
the experimental data with theoretical predictions to get some insight into the pinning 
mechanism in MgB 2 . Using A oc (1 — t 4 ) -1 / 2 and £ oc [(1 + t 2 )/(l — t 2 )] 1 / 2 , Griessen et al. 
have found that for ^/-pinning the critical current density in the single vortex pinning regime 
j sv oc (1 - t 2 fl 2 {\ + t 2 )" 1 / 2 , while for <5T c -pinning j sv oc (1 - t 2 ) 7 / 6 (l + 1 2 ) 5 / 6 . Inserting all 
these expressions into Eq.(||), we have 



is obvious. 




(5) 



for £T c -pinning, and 




(6) 
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for <5Z-pinnirig. The lines corresponding to Eqs.fj^) and (||) are indicated as £T c -pinning and 
^/-pinning respectively in Fig. |^. The central result of this Letter is the remarkably good 
agreement found between 5 s b and the corresponding <5T c -pinning line in the figure. In sharp 
contrast, the <5Z-pinning line shown in Fig. |^ is in total disagreement with the experimental 
data. 

Having derived the crossover fields B s ^ and _B t h, we now reconstruct the B-T phase 
diagram shown in Fig. |||. The final B-T phase diagram is shown in Fig. || The vortex solid 
region is divided into three smaller regions. Dingle vortex pinning governs the region below 
_B s b, between B s ^ and B t h, small bundle pinning becomes dominant, while between B t h and 
.6^,., thermal fluctuations are more important. Large flux bundle pinning is not observed in 
MgB 2 , but may be concealed by the thermal fluctuation effects. 

In summary, we have found strong evidence for £T c -pinning, i.e. pinning via the spa- 
tial fluctuations in the transition temperature, in the new superconductor MgB2, while 
^/-pinning i.e., pinning via the spatial fluctuations of the charge carrier mean free path, is 
not observed. The B-T phase diagram of the MgB 2 sample has been derived, showing that 
at low fields below B s ^ the system is dominated by single vortex pinning and changes to 
smaller bundle pinning when B > B s ^. When B > B t ^, this region in the vortex solid area 
is dominated by thermal fluctuations. The irreversibility line may result from the giant flux 
creep effect. 

The authors would like to thank the Australian Research Council for financial support. 
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FIGURES 

FIG. 1. Magnetization hysteresis loops of the MgB2 bulk sample taken every 2 K in the 14-36 
K range. Results at lower temperatures are not shown because of large flux jumping. Inset 
shows the magnetization hysteresis loop of a pressed MgB2 sample at 5 K, showing large reversible 
magnetization from the surface current. 

FIG. 2. Critical current density j c calculated from the Bean critical state model, indicated by 
different symbols for different temperatures. Solid lines are fitting curves using Eq.(|j). 

FIG. 3. B-T phase diagram of the new superconductor MgB2- Open circles represent the 
irreversibility line obtained from Fig. [2l and the solid line is a fit to B irr = 5.2(1 — t 2 ) 3 / 2 . Solid 
circles represent the B C 2{T) line from Takano's data using resistive measurements. The dashed line 
is just a guide to the eye. 

FIG. 4. Critical current density at 24 K in a double-logarithmic plot of — log[j/j(-B = 0)] 
versus the applied field. The solid line is a fit using Eq.([|). B s ^ indicates the crossover field 
from single vortex pinning to small bundle pinning, while B t ^ is the crossover field to the thermal 
fluctuations dominated regime. Inset shows the j c (B) curve of the pressed MgB2 sample with a 
large contribution from the surface current, which shows a B~ 12 behavior. 

FIG. 5. Temperature dependence of the crossover field B s ^. The <5T c -pinning line corresponds 
to Eq- (|5j) , which is in agreement with the experimental data, while The ^-pinning line corresponds 
to Eq.(^), which is not in agreement with the experimental data. 

FIG. 6. B-T phase diagram of the new superconductor MgB2- B s b(T) is the fitting curve 
of Eq.(^) to the experimental data (see Fig. ||). B[ TT (T) is the fitting curve of Eq.([l]) to the 
experimental data (see Fig. ||). B^ is the crossover field to thermal dominant region (see Fig. |5|). 
Again the B C 2(T) line is taken from Takano's data using resistive measurements. 



S 



REFERENCES 



[1] J. Nagamatsu, N. Nakagawa, T. Maranaka, Y. Zenitani, and J. Akimitsu, Nature 410, 
63 (2001). 

[2] D. C. Larbalestier, M. O. Rikel, L. D. Cooley, A. A. Polyanskii, J. Y. Jiang, S. Patnaik, 
X. Y. Cai, D. M. Feldmann, A. Gurevich, A. A. Squitieri, M. T. Naus, C. B. Eom, E. 
E. Hellstrom, R. J. Cava, K. A. Regan, N. Rogado, M. A. Hayward, T. He, J. S. Slusky, 
R Khalifah, K. Inumaru, and M. Haas, Nature 410, 186(2001). 

[3] Y. Bugoslavsky, G. K. Perkins, X. Qi, L. F. Cohen, and A. D. Caplin, Nature 410, 563 
(2001). 

[4] D. G. Hinks, H. Claus, and J. D. Jorgensen, Nature 411, 457 (2001). 

[5] Y. Bugoslavsky, L. F. Cohen, G. K. Perkins, M. Polichetti, T. J. Tate, R. Gwilliam, 
and A. D. Caplin, Nature 411, 561 (2001). 

[6] S. Jin, H. Mavoori, C. Bower, and R. B. van Dover, Nature 411, 563 (2001). 

[7] C. B. Eom, M. K. Lee, J. H. Choi, L. Belenky, X. Song, L. D. Cooley, M. T. Naus, 
S. Patnaik, J. Jiang, M. Rikel, A. Polyanskii, A. Gurevich, X. Y. Cai, S. D. Bu, S. E. 
Babcock, E. E. Hellstrom, D. C. Larbalestier, N. Rogado, K. A. Regan, M. A. Hayward, 
T. He, J. S. Slusky, K. Inumaru, M. K. Haas, and R. J. Cava, Nature 411, 558 (2001). 

[8] J. S. Slusky, N. Rogado, K. A. Regan, M. A. Hayward, P. Khalifah, T. He, K. Inumaru, 
S. Loureiro, M. K. Haas, H. W. Zandbergen, and R. J. Cava, Nature 410, 343 (2001). 

[9] P. C. Canfield, D. K. Finnemore, S. L. Bud'ko, J. E. Ostenson, G. Lapertot, C. E. 
Cunningham, and C. Petrovic, Phys. Rev. Lett. 86, 2423 (2001). 

[10] D. K. Finnemore, J. E. Ostenson, S. L. Bud'ko, G. Lapertot, and P. C. Canfield, Phys. 
Rev. Lett. 86, 2420 (2001). 

[11] S. L. Bud'ko, C. Petrovic, G. Lapertot, C. E. Cunningham, and P. C. Canfield, M. H. 

9 



Jung, A. H. Lacerda, Phys. Rev. B 63, 220503 (2001). 

[12] Herbert Schmidt, J. F. Zasadzinski, K. E. Gray, and D. G. Hinks, Phys. Rev. B 63, 
220504 (2001). 

[13] Y. Takano, H. Takeya, H. Fujii, H. Kumakura, T. Hatano, K. Togano, H. Kito, H. Ihara, 
|cond-mat/0l02l67| . 



[14] H. H. Wen, S. L. Li, Z. W. Zhao, Y. M. Ni, Z. A. Ren, G. C. Che, Z. X. Zhao, |cohd 



mat/0103521 . 



[15] Amish G. Joshi, C. G. S. Pillai, P. Raj, S. K. Malik, |cond-mat/0 103302 



[16] R. Griessen, Wen Hai-hu, A. J. J. van Dalen, B. Dam, J. Rector, and H. G. Schnack, 
Phy. Rev. Lett. 72, 1910 (1994). 

[17] G. Blatter, M. V. Feigel'man, V. B. Geshkenbein, A. I. Larkin, V. M. Vinokur, Rev. 
Mod. Phys. 66, 1125 (1994). 

[18] Roger Wordenweber, Rep. Prog. phys. 62, 187 (1999). 

[19] E. V. Thuneberg, Cryogenics 29, 236 (1989). 

[20] E. V. Thuneberg, J. Kurkijarvi, D. Rainer, Phys. Rev. Lett. 48, 1853 (1982). 
[21] C. J. van der Beek, and P. H. Kes, Phys. Rev. B 43, 13032 (1991). 
[22] P. H. Kes, Physica C 185-189, 288 (1991). 

[23] S. X. Dou, X. L. Wang, J. Horvat, D. Milliken, E. W. Collings, and M. D. Sumption, 
|cond-mat/010232q . 

[24] M. J. Qin, X. L. Wang, S. Soltanian, A. H. Li, H. K. Liu, and S. X. Dou, Phy. Rev. B 
64, 060505 (R) (2001). 

[25] K. A. Miiller, M. Takashige, and J. G. Bednorz, Phys. Rev. Lett. 58, 1143 (1987). 



10 



[26] Y. Yeshurun, and A. P. Malozemoff, Phys. Rev. Lett. 60, 2202 (1988). 
[27] M. Tinkham, Phys. Rev. Lett. 61, 1658 (1988). 
[28] E. H. Brandt, Rep. Prog. Phys. 58, 1465 (1995). 

[29] A. I. Larkin, and Yu. N. Ovchinnikov, J. Low Temp. Phys. 34, 409 (1979). 



11 




Fig. 1 by Qin et al. 
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Fig. 3 by Qin et al. 
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